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ABSTRACT 


The potential of a suitably designed and economically manufactured 
Stirling engine as the energy conversion subsystem of a paraboloidal 
dish-Stir ling solar thermal power module has been estimated. 

Results obtained by elementary cycle analyses have been shown to match 
quite well the performance characteristics of an advanced kinematic Stirling 
engine, the United Stirling P-40, as established by current prototypes of the 
engine and by a more sophisticated analytic model of its advanced derivative. 

In addition to performance, brief consideration has been given to other 
Stirling engine criteria such as durability, reliability, and serviceability. 
Production costs have not been considered here. 


in 



FOREWORD 


A few years ago, an investigation was undertaken as part of a comparative 
study of promising advanced power conversion subsystems suitable for use in 
paraboloidal dish solar thermal power systems. The results of this investigation 
are reported herein. These results have been recently reviewed and are 
considered currently applicable. Ultimately, the Stirling engine and other 
selected types of engines will be compared on the basis of performance, 
durability, reliability, serviceability, and cost criteria. The potential 
rewards and risks associated with future technological and economic developments 
that affect these criteria must also be reassessed in each case. 
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SECTION I 


INTRODUCTION 


A few years ago, an intensive effort was made to determine the potential 
of a suitably designed and economically manufactured Stirling-cycle engine as 
the energy conversion subsystem of a point-focusing, cavity-receiver solar 
thermal power module. Currently, the module of greatest interest comprises a 
paraboloidal collector that is pointed continuously at the sun and carries a 
cavity receiver and power conversion unit at its focus; however, the 
applicability of the Stirling-cycle engine is not limited to this particular 
configuration. 

As with the Carnot cycle, the ideal Stirling cycle is reversible and, 
thus, has the highest thermal efficiency attainable when operated between 
prescribed temperatures of source and sink. Yet, as is also true for the 
Carnot cycle, no practical engine has yet been devised that operates on the 
ideal Stirling cycle. 

Some practical engines operate on one or the other of a variety of 
cycles that approximate the ideal Stirling cycle. Unfortunately, each of 
these practical engines incorporates apparently unavoidable design features 
that limit the compression ratio of the engine and introduce thermodynamic 
irreversibilities and mechanical losses; as a result, the thermal efficiency 
of these practical engines is universally less — frequently much less — than 
that of a reversible engine operating between the same temperatures of source 
and s ink . 

An equally troublesome set of limitations that has persistently retarded 
the development of a practical Stirling-cycle engine consists of (1) the 
sliding seals required for pistons and the sliding or rolling seals required 
for piston rods, (2) the oil contamination that fouls regenerators in 
particular and the entire working space of the engine in general when piston 
and/or piston rod seals leak, (3) the high cost of fabricating regenerators 
and heaters that are highly effective at the moderately high cyclic speeds 
required to obtain a high specific power-to-volume or power-to-weight ratio, 
(4) the high mean working pressure required for a compact engine, (5) the 
kinds and quantities of metals and/or ceramics required to allow operation at 
source temperatures high enough to provide an engine power density and 
efficiency competitive with internal combustion engines, and (6) the 
relatively large cooling system required to reject the heat from this 
closed-cycle engine as compared with open-cycle engines. 

Certain attractive features of a practical engine operating on a cycle 
that approximates the Stirling cycle have served to stimulate a sustained 
interest in this concept in the face of all the problems noted above. Some of 
these features are (1) potentially high thermal efficiency, (2) quietness and 
smoothness of operation, and (3) an ability to utilize a variety of fuels with 
relatively less resultant atmospheric pollution than that generated by a 
comparable internal combustion engine. 
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The following sections of this report discuss (1) the Stirling and 
quasi-Stirling cycles and simulation of these cycles, (2) the limitations that 
currently circumscribe the practicality of so-called Stirling engines and the 
probability that these limitations may be relaxed or removed, and (3) an 
attempt to assess the potential of the practical quasi-Stirling engine as the 
power conversion subsystem for a point-focusing cavity-receiver, solar thermal 
power module. 


A. OBJECTIVES 

In view of the remarks presented in preceding paragraphs, the following 
objectives were defined for this investigation: 

(1) Determination of the principal parameters affecting Stirling 
engine performance and assessment of the sensitivity of 
performance to variations in these parameters. 

(2) Estimation of the potential of the kinematic Stirling engine in 

terms of the following criteria: performance, durability, 

reliability, and serviceability. Mass production costs have been 
addressed in a separate study (Ref. 1). (Thermodynamic cycle 
analysis, examination of the above-named criteria as they apply to 
recent and currently existing Stirling engines, and projection of 
these criteria for future engines that incorporate promising new 
theoretical and technological innovations have each contributed to 
the results presented here.) 

(3) Suggestion of a development plan that offers the probability of 
considerable improvement in kinematic Stirling engines as measured 
by the stated criteria. 

B . BACKGROUND 

Earlier studies (Refs. 2, 3, and 4) noted the potentially high 
performance capability of Stirling engines; however, these studies also noted 
the relatively undeveloped state of this type of engine as compared, for 
example, with those of the Rankine- and Brayton-cycle engines. Of the two 
general types of Stirling engine, kinematic and free piston, the free piston 
type was assumed to have the higher ultimate potential, but to be less 
well developed than the kinematic type. Thus, the kinematic type Stirling 
engine was chosen as the more viable power conversion unit for initial 
paraboloidal dish-Stir ling solar thermal power plants (Ref. 5). 

Two existing kinematic Stirling engines, viz., the United Stirling P-40 
and the Phillips 1-98, were adapted to serve as the power conversion sub- 
systems in first-generation paraboloidal dish-Stir ling solar thermal power 
module designs. Even these engines could not be considered to be fully 
developed at the time this study was done; thus, they served in the study 
reported here only as the base from which the criteria for future engines have 
been projected. 
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SECTION II 


THE STIRLING CYCLE — DISCUSSION AND ANALYSIS 


A. DESCRIPTIONS AND DEFINITIONS 

Two rather different general configurations for Stirling-cycle engines 
have been developed. The kinematic configuration, like the usual Otto- or 
Diesel-cycle engine, converts the linearly oscillatory motion of the power 
piston to steady rotational motion of a shaft. In the free-piston configura- 
tion, the oscillatory linear motion of the power piston can be directly 
imposed on a coaxial extension of the piston rod that, in turn, drives some 
linearly oscillating load, e.g., the armature of a linear electric alternator. 

The kinematic configuration is represented by two major versions, viz., 
the dual-piston version and the piston-displacer version. In turn, two 
variants of the piston-displacer version exist. Figure 2-1, taken from 
Reference 6 , defines these three kinematic configurations. 

Figure 2-2 presents the pressure-volume diagrams for the simple 
regenerative Stirling cycle and for the regenerative Stirling cycle with 
regenerator dead volume. Point 4' of the latter cycle has been made coinci- 
dent with point 4 of the former in order to provide a definite basis for 
comparison in the following discussion. 

Figure 2-3 presents the sequence of processes constituting the cycle for 
a dual-piston, double-acting (Rinia-type) Stirling engine. Such an engine 
usually comprises four cylinders connected in a ring arrangement. 

Figure 2-4 presents the sequence of processes constituting the cycle for 
a Beta-type piston-displacer Stirling engine. 

Elementary thermodynamic analysis of the ideal Stirling cycle is 
relatively simple; however, because the ideal cycle only crudely represents 
the thermodynamic cycle on which any practical Stirling engine operates, the 
results of the elementary analysis are of little practical significance. In 
order to provide an introduction to the development of a second, more 
elaborate model and to provide a reference with which that second model can be 
compared, the indicated efficiency and work performance parameter of the 
simple regenerative Stirling cycle is presented here. 

The ideal Stirling cycle consists of a series of four discrete processes 
(cf. Figure 2-2): an isothermal compression from thermodynamic state (p 3 , 

Vf, T 3 ) to thermodynamic state (p 2 > T 2 = Tj), a constant-volume 

heat addition from thermodynamic state (p 2 > V 3 > T 2 ) to thermodynamic state 
(p 3 , ^3 = ^2> T 3 ), an isothermal expansion from thermodynamic state (p 3 , V 3 , T 3 ) 
to thermodynamic state (P 4 , V 4 , T 4 = T 3 ) , and a constant volume heat rejection 
from thermodynamic state (P 4 , V 4 , T 4 ) back to the original thermodynamic state 
<P1, V X = V 4 , Tjl). 
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DUAL-PISTON VERSION 


PISTON-DISPLACER VERSIONS 



ALPHA-TYPE 


LEGEND: 

H = HEATER 
R = REGENERATOR 
C = COOLER 

1 = EXPANSION SPACE 

2 = COMPRESSION SPACE 



BETA-TYPE 



Figure 2-1. Stirling Engine Kinematic Configurations 






Figure 2-2. Simple, Perfectly Regenerative Stirling Cycle: 1—2— 3— 4—1 ; 

Perfectly or Imperfectly Regenerative Stirling Cycle with 
Regenerator Dead Volume: 1 ' —2 * —3 1 —4 * — 1 * (Note: V 4 ' = V 4 
and P 4 ' = P 4 .) 


* 


When regeneration is perfect, the thermodynamic efficiency of the ideal 
Stirling cycle is equal to that of the Carnot cycle when both cycles are 
operated between the same temperatures of source and sink. This efficiency is 
expressed by the familiar formula, 


n s (ideal) ^Carnot T 

ri 


(1) 


This is shown in the discussion of the simple regenerative Stirling cycle that 
appears in the following section of this report. 

The ideal Stirling cycle analysis presumes the ability to satisfy 
several thermodynamic conditions that are not physically realistic. Among 
these are the following: 

(1) During the expansion or compression process, all of the 

thermodynamic medium is assumed to be at temperature Tjj or Tq, 
respectively. Thus, the ideal cycle fails to allow for those 
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volumes of the dead space within the cooler at temperature Tc 
during the expansion process, within the heater at temperature 
Tfl during the compression process, or within the regenerator at 
varying temperatures intermediate to Tjj and Tq during either 
process . 

Furthermore, unless elaborate efforts are made to maintain a more 
or less uniform, constant temperature of the medium that is 
undergoing expansion or compression within the hot and cold 
spaces, the gas in these spaces will be subject to processes that 
are more nearly adiabatic than isothermal. 

(2) During the constant-volume heat addition and removal processes of 
the cycle, all of the thermodynamic medium is assumed to be at one 
common temperature at any instant, even though some of the working 
fluid must again reside in the cooler at temperature Tc, in the 
heater at temperature Tg, and in the regenerator at various 
intermediate temperatures. 

(3) No actual materials are perfect insulators or, at ordinary 
temperatures, perfect conductors. Because any of the usual 
thermodynamic media must be contained at elevated temperatures and 
pressures during at least part of any practical thermodynamic 
cycle, boundaries of finite strength and corresponding thickness 
are required. Heat leaks and heat transfers across finite 
temperature differences occur at these boundaries; these processes 
are thermodynamically irreversible and therefore result in reduced 
cycle efficiency. 

(4) All flow processes within the working space are assumed to occur 
without loss, i.e., the thermal conductivity of the fluid is 
assumed to be infinite, yet the closely related viscosity of the 
medium is assumed to be zero. Unless the thermodynamic processes 
of the cycle all proceed at infinitesimal rates, the ideal cycle 
cannot be realized when any real thermodynamic medium is used. 


B. ANALYSIS 

1. Evaluation of the Simple Regenerative Stirling Cycle 

a. Zero Dead Volume; Perfect or Imperfect Regeneration . The 
work accomplished by a reversible, isothermal process is given by the formula, 


* 


W = mRT In 


1 


( 2 ) 
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Applied to the expansion and compression processes of the Stirling 


cycle , 


V V 

4 1 

W = mRT In — = mRT„ In 

t* rl V3 n v 2 


( 3 ) 


and 


4 


V V 

W c = ®RT C In ^ - mRT c In ^ , 


( 4 ) 


respectively. The net work per cycle is, thus, 



( 5 ) 


The thermal energy furnished to the engine per cycle consists of the 
thermal equivalent of the work of expansion plus the constant-volume heat 
addition required to compensate for imperfect regeneration: 


Qi n - mC v (T H - T R ) + W E . (6) 

Defining regenerator effectiveness as 



After algebraic manipulation, Q^ n can now be expressed as: 



( 8 ) 
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Cycle efficiency can now be determined: 


n 


w 


net 


in 


or 



( 9 ) 


Obviously, when regeneration is perfect, E = 1 and 


n 



n 

Carnot 


( 10 ) 


Defining the cycle work performance parameter as 


P 

w 


W 


net 


n (v i - V 


on 


and substituting appropriately, 

I 

( 12 ) 

Note that here the minimum cycle pressure, pj, has been employed in the 
definition of P w . 
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b . Finite Dead Volume; Perfect or Imperfect Regeneration . In a 
practical Stirling-cycle engine, that part of the volume occupied by the 
working medium but not swept by pistons is termed the "dead space." Dead 
space comprises the piston clearance volumes, the heater, regenerator, and 
cooler net internal volumes, and the volume of the connecting ducts. A thermo- 
dynamic analysis that is slightly more elaborate than the most elementary 
retains the concept of a series of distinct thermodynamic processes but allows 
consideration of the dead space as well as of regenerator effectiveness. This 
second method of analysis provides a considerably more accurate model for the 
performance of a practical Stirling engine. The trends of indicated engine 
efficiency and indicated work performance parameter as functions of regenerator 
effectiveness and/or dead volume are correctly indicated. Still, the analysis 
employs lumped parameters and distinct thermodynamic processes, so, at best, 
the resulting model must be calibrated to the corresponding practical Stirling 
engine before useful quantitative results can be obtained. 

This second model has been the primary analytic tool of this study. The 
much greater complexity of the next significantly higher level of analysis, 
i.e., the second-order methods, prevented its application to the problem at 
hand at this time; nevertheless, the glaring weaknesses of the model employed 
here must be admitted, and the results obtained from it must be circumspectly 
evaluated . 


2. Regenerative Stirling Cycle with Regenerator Dead Space 

First, consider the discharge of the regenerator. Assume that the 
regenerator discharge process results in the advance from one end of the 
regenerator to the other of a discontinuous drop in temperature from the cycle 
source temperature to the cycle sink temperature as shown in Figure 2-5. 
Further assume that a constant average pressure exists in the regenerator 
during this process (this as an approximation). 

Under these approximations, the time-averaged mass of gas contained in 
the regenerator void volume resulting from this discharge process will be 


» 



(13) 


where T, the instantaneous average temperature of the gas in the regenerator 
void volume, is given by 


T 


T C + 




(14) 
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The charging of a regenerator can be considered in an exactly 
corresponding way and results in the same expression for Tr. 

Several subtle but appreciable approximations are hidden in this 
derivation of a representative time-averaged regenerator temperature; some of 
these are 

(1) It is assumed that the regenerator discharging or charging process 
is complete. This cannot be precisely true. 

(2) Use of this average regenerator temperature implies that the 
regenerator is always either charging or discharging. This is 
untrue because, at certain phases of the cycle, the thermal energy 
in the regenerator remains virtually constant at either the 
extreme high level (fully charged) or the extreme low level 
(completely discharged). 

(3) The thermal capacity of either the heating or cooling gas and the 
heat transfer coefficient between either heating or cooling gas 
and the regenerator matrix are all finite. 

(4) The step increase or decrease in temperature in the regenerator is 
dispersed by axial conduction in the regenerator matrix. 

(5) Part of the thermal energy stored in the regenerator as a result 
of the charging process is lost by thermal leakage . 

(6) A constant average pressure during the charging or discharging was 
assumed in obtaining the log-mean temperature difference for Tr; 
however, in the following development, Tr will be employed in a 
formula for instantaneous pressure during the same charging or 
discharging process. The effect of this approximation is unclear. 

A second, more conventional approach to determining the time-averaged 
mass and mean temperature of the gas in the regenerator utilizes the analysis 
of H. Hausen (Ref. 7). A translation of the most pertinent parts of this 
reference is appended. Hausen showed that the mean spatial distribution of 
the temperature in the matrix is linear and constant under the following 
assumptions : 

(1) Equal thermal capacities and half-cycle times for the heating and 
cooling gases. 

(2) Negligible end effects (end effects are small in the usual 
regenerator) . 

(3) Large thermal capacity of the matrix relative to those of the 
gases . 

(4) Negligible axial conduction in the matrix. 
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(5) Constant heat transfer coefficient between the gases and the 
matrix . 

Under these conditions, Hausen also showed the temporal mean temperature 
difference between the heating gas and the matrix to be (1) equal to that between 
the matrix and the cooling gas and (2) constant. 

If the further assumption is made that the heat transfer coefficient 
between either gas and the matrix is very large, the temperature difference 
between gas and matrix will be relatively small. The temporally constant 
distribution of the gas mean temperature can then be written as 


T 




(18) 


Assuming a constant average pressure in the regenerator, as was done in 
the first case, it can again be shown that 


m R 


P V 


R 


R T. 


R 




where 


T 


R 



This same result has been presented previously by Martini (Ref. 6) and 
Berchowitz (Ref. 8). Comparison of these two analyses based on radically 
different assumptions shows the results for Tp to be identical and strongly 
reinforces the author's belief that Tr is probably relatively insensitive to 
the choice of regenerator model. 

In spite of the inaccuracies imposed by using a constant log-mean 
temperature difference to represent a time-dependent average regenerator 
temperature, the indicated efficiency and work performance parameter of the 
resulting model have been shown a posteriori to be reasonably accurate when 
the limited parameters of the model are determined by calibration of the model 
relative to experimental results and other more elaborate models . 
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During the isothermal expansion process in this more complicated model, 
the working gas inventory is represented by 


R \ T H V 


(Note that no allowance has been made in this model for any void volume at the 
sink temperature, even though practical Stirling engines require some 
clearance volume at the cold end of the working space and, at least in 
contemporary engines, some volume in the cooler.) 

Solving for p', 


Determining the work of expansion 


/ V 1 

V 

^ + J 

T„ T_ 


m'RT In 
H. 


V V 

JL + _5 

t h t r 

V V 
- 2 + 

t h t r 


m'RT In 
H 


T V 

T V 
H 1 
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Determining the work of compression in a similar way, 



TV TV 
-1 12l + — — 
TV TV 
H 1 H 1 

T TV 
3 . _C _R 

T TV 
H H 1 


Determining the thermal energy added per cycle, 


Q'in =m ' C v (T H-V +W E 



( 22 ) 


( 23 ) 

1 
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The indicated efficiency of this model can now be determined 



( 24 ) 


From the above equation, it can be seen that if 


E = 1 and 



> 


then 


n' = 1 - 


= n 


H 


Carnot 


Also it can be seen that when 



> 


2-15 


then 


= "Carnot ‘ 0 ' 


The latter case is, of course, trivial. 

In all other practical cases, even when E = 1, if 



< 1 , 


and 


then 



< 1 


1 



0 , 


n' 


< ^Carnot * 


The work performance parameter for the regenerative Stirling cycle with 
regenerator dead space is defined as 


P' 

w 


W' 


net 


P ; <v t 


V 


(25) 


Thus, the form is exactly the same as that employed in the case of the simple 
regenerative Stirling cycle. 
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As in Reference 6 , the pressure, P 4 ', at the end of the isothermal 
expansion stroke in the regenerative Stirling cycle with regenerator dead 
space has been made equal to the pressure, P 4 , of the same point in the 
simple regenerative Stirling cycle. The rationale for this choice is that 
P 4 and p 4 ' closely approximate the temporal mean pressure of the 
corresponding model, and, for a specified Tj{, creep strength based on 
temporal mean pressure is the criterion that determines the wall thicknesses 
of the cylinder head, heater tubes, and other hot-end pressure-retaining 
components . 

It is possible to show by a rather tedious analysis that, given 
P4' = P4> 


and 



m 


1 + 


V /V 
R /V l 

T /T 
R' H 


> 




T TV 

T TV 

H H V 1 





T V 

J* 

T + V 
H 1 


TV TV 
3 _2 _£ _R 

T V + T V 
H 1 H 1 




H 


T V 
_R _2 

T V 
H 1 



(26) 


(27) 


(28) 


(29) 
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Whether regeneration is perfect or imperfect, the ratio of the net works 
performed, henceforth denoted simply by W'/W, is given by 



In 


H 

T V 
_R _2 

‘T V 
H V 1 



_ V TV 
_2 + _C _R 

‘TV + T V 
H 1 H 1 


(RR 


( 30 ) 


v 


Now, resuming the discussion of P w ', comparison shows that 


P 1 

w _ W 1 /W 
P w ' p l' /p l 


(31) 


and , 


substituting for W'/W and pi'/pi 


leads to the relation 


P ' 
w 


P 

w 


1 + 



R 


In 


H 


T V 
_R _2 

■ T V 
H 1 


T T. V„ 

R | C R 

T TV 
C / H H V 1 

T in l T V TV 
H V R 2 . C R 


V, 


■ T V 
H 1 


T V 
H 1 



(32) 


The maximum working volume of the regenerative Stirling engine with 
regenerator dead space is V max = V^ + Vr, while the minimum working 
volume of the engine is V m £ n = V 2 + Vr. It follows from these two 
equations that 


V R 

V R 

V 1 

or 

1 < 
11 

V R 

V 

max 

1 

V 

V, 

V„ 

max . R 

V 1 


1 

L - R 
V 

max 


( 33 ) 
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and 


V 

max 

V . 
mm 


1 + 




or 


V. 


1 - 


min 


max 


min 


mm 


(34) 


Of course, still other volume ratios could be defined; however, the 
important point to note here is that definition of the simple regenerative 
Stirling cycle requires the specification of k, Tq/T r , X^/V^, Vp/V^, and E. 

In the case of the simple cycle, V 2 /V 1 = V ra £ n /V max ; in the case of the more 
complicated cycle, the two volume ratios V 2 /V 1 and V R /V£ can be replaced by 
^max^min anc * ^R/^min ^ desired. 

The effect of regenerator dead space on cycle performance is perhaps best 
illustrated by a specific example. If it is assumed that k = 1.4, T^/T^ = 1/3, 
V 2 /V 1 = 1/2, and E - 1, the results presented in Figure 2-6 are obtained. 

In the preparation of this graph, the assumption of equal average 
working space pressure at the end of the expansion strokes of the engines 
without and with regenerator dead space was made. Note that the work ratio, 

W'/W, and the efficiency ratio, n' / n c arn ot » decrease nearly linearly with 
increasing ratio of regenerator volume to maximum working space volume. 

The effects of variations in k, T(j/T H , V max /V m £ n , V R /V m £ n , and E 
on the regenerative Stirling cycle with regenerator dead space have been explored 
parametrically. Results of this parametric evaluation of the model are presented 
in Figures 2-7 through 2-11. In Figure 2-11, represents the mechanical 
efficiency of the engine. 

Because of the heat transfer and thermal energy storage requirements of 
the heater, regenerator, and cooler components coupled with the requirement 
for reasonably low pressure losses in these same components, there exists a 
practical maximum limit to volume compression ratio, Vmav/Vnn' p . Walker 
(Ref. 9) has estimated this limit to be V max /V ra £ n =: 2.5. The corresponding 
maximum for the ratio of regenerator dead volume to minimum volume is 1.0; 
however, this can occur only if regenerator dead volume and minimum volume are 
equal! This is not a practical possibility because of (1) the necessity for 
displacing the minimum volume of working gas through the regenerator in order 
to discharge the regenerator and, thus, to heat the gas at constant volume 
during this phase of the thermodynamic cycle and (2) the required finite 
volumes of the heater and the cooler. After considering these practical 
requirements, it was concluded that a practical maximum value for the ratio of 
regenerator dead space to minimum volume in a conventional Stirling engine is 
approximately one-half. 
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Figure 2-6. Effect of Dead Volume on Discrete-Process Stirling Cycle 
Performance 
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Figure 2-7. Normalized Work Performance Parameter as a Function 
of Dead Space to Minimum Working Space Ratio for a 
Discrete-Process Stirling Cycle 
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Figure 2-9. Indicated Efficiency as a Function of Maximum to Minimum 
Working Space Ratio for a Discrete-Process Stirling Cycle 



INDICATED EFFICIENCY (v) 


REGENERATOR EFFECTIVENESS 


t r- t c 


t h- t c 


V /V 
dead' min 


k = 1.4 

k = 1.67 

W 3 - 5 

V /V . =2.5 

max mm 


Figure 2-10. Effect of Regenerator Effectiveness on Stirling Cycle 
Efficiency for a Discrete-Process Stirling Cycle 
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EVALUATION OF RESULTS OBTAINED FROM THE MODEL 


The pronounced effect of dead volume to minimum volume ratio, 

^dead/ v min> on work performance is indicated in Figure 2-7. At 
^dead/^min = 0*5, an almost one percent degradation in normalized work 
performance parameter, P w '/P w > results from each one percent increase in 
^dead/^min* This degradation is more pronounced at smaller values of 
^dead/^min anc * smaller values of the maximum to minimum volume ratio, 

V max /V m i n , although, over the ranges of these parameters shown in the 
figures, the effect due to V max /V m £ n is obviously less than the effect due 
to V dead/ V min* 

As shown in Figure 2-8, the sensitivity of indicated efficiency, n , to 
^dead/^min on ly about two-thirds as large as the sensitivity of P w '/P w 
to ^dead/^min i- n the vicinity of the representative engine points; however, the 
marked increase in this sensitivity at increasing values of should be 

noted . 


The effect of maximum to minimum working space ratio, V max /V m; £ n , 
on n is presented in Figure 2-9; the sensitivity of n to V max /V m: [ n is not 
large in the vicinity of the engine points used in the example, but this 
sensitivity is greater at smaller values of V max /V m ^ n> 

Because of the scales selected for the ordinate and abscissa of 
Figure 2-10, the sensitivity of n to regenerator effectiveness, E, may not 
appear impressive; however, in the vicinity of the sample engine points, a 
degradation of more than one percent in n results from a one percent decrease 
in E when the working medium is a diatomic gas and a two-thirds percent 
decrease in n with a one percent decrease in E when the medium is a monatomic 
gas . 


From Figures 2-8, 2-9, 2-10, and 2-11, the observant reader may have 
noted that the indicated engine efficiency obtained by using a monatomic 
medium is apparently better than that obtained by using a diatomic medium when 
regeneration is less than perfect, i.e., E<1. Actual Stirling engines almost 
universally employ either helium (monatomic) or hydrogen (diatomic) as the 
working medium, and, in the same engine, hydrogen often results in a slightly 
higher efficiency than does helium. Both hydrogen and helium are employed 
because of their low molecular weights and associated excellent heat transfer 
and low flow loss characteristics; in these respects, hydrogen is superior to 
helium. However, the analytic model employed herein does not account for heat 
transfer and flow loss effects, and, thus, helium (being monatomic) exhibits 
its thermodynamic advantage over hydrogen when regenerator effectiveness is 
less than one. There appears to be a consensus that, when all factors are 
considered, either of these two gases in a suitably optimized variant of the 
same engine would result in about the same efficiency as the other. 

Engine efficiency (indicated or brake, as noted on particular curves) is 
presented as a function of working gas high temperature in Figure 2-11. 
Regenerator effectiveness, E, and the ratio of specific heats of the working 
medium, k, also appear as parameters in this figure. Note that this figure 
pertains to one sample engine (V max /V m i n = 2.5 and V c j ea< j/V m -[ n = 0.5) 
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that is rejecting heat to a thermodynamic sink at 140°F. The principal 
characteristics of the Stirling engine that are indicated by Figure 2-11 are 
the following: 

(1) High regenerator efficiency becomes increasingly important as the 
working gas high temperature is increased . 

(2) When regenerator effectiveness is imperfect, the thermodynamic 
advantage of helium over hydrogen that is indicated by the model 
employed here becomes increasingly important as the working gas 
high temperature is increased. 

(3) The magnitude of the kinematic Stirling engine brake efficiency 
predicted for a regenerator effectiveness of 0.95 at 1500°F 
working gas high temperature is in close agreement with the 
performance predicted for the improved United Stirling P-40 and 
P-75 engines, viz., n~40%. Furthermore, insofar as can be 
determined at present, the trend of this efficiency with 
increasing working gas high temperature approximates the predicted 
performance of advanced technology derivatives of these same 
engines . 


D. HIGHER-ORDER METHODS OF STIRLING- CYCLE ANALYSIS 

In order to consider the more realistic quasi-Stirling cycle 
corresponding to the sinusoidal motions of the pistons and the resulting 
overlapping thermodynamic processes, so-called first-order and second-order 
methods of analysis have been developed, cf. Ref. 6. These methods may be 
either numerical or analytic; however, numerical analysis allows the more 
general consideration of regenerator storage characteristics. As customarily 
applied, only the work performed per cycle and the work performance parameter 
are calculated on a strictly analytic basis. The indicated cycle efficiency 
is estimated by applying what is termed an experience factor to the Carnot 
efficiency corresponding to the given temperatures of source and sink, and the 
cycle brake efficiency is estimated by applying a second experience factor, 
the mechanical efficiency, to the indicated cycle efficiency. The two 
experience factors vary with the type of Stirling engine, its method of 
construction, the operating temperatures, etc.; thus, these factors are 
heavily dependent on experience and must be determined from existing data. 

While this method has been applied in the course of the current study, 
no results are presented here because of the tenuously empirical nature of the 
method . 

Also in this study, an initial attempt was made to employ the very 
sophisticated third-order numerical program published by Tew, et al., 

(Ref. 10) for analysis of the Stirling engine; however, this attempt was 
terminated for the following reasons: 

(1) As written, the program was applicable only to piston-displacer- 
type engines . 
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( 2 ) 


Numerical values for many of the engine design parameters were 
inserted directly into the program rather than through input 
lists, and the values inserted were uniquely representative of the 
GM GPU- 3- 2. 

(3) As a result of two keypunch errors at JPL that could not be 

detected by the computer diagnostic program, correct operation of 
the Stirling program was obtained only after tedious and time- 
consuming examination and too late to be utilized in this study. 

Of the models described above, only the model of Tew, et al., and other 
similar sophisticated numerical models can emulate Stirling engine performance 
and efficiency without a heavy dependence upon empirically established 
correction factors. 


E. FURTHER REMARKS 

No practical thermodynamic engine has ever demonstrated an efficiency 
equal to that of its corresponding idealized cycle; in fact, the practical 
utility of thermodynamic engines has, in the past, depended as much on the 
ability of the actual engine to approach the corresponding ideal as on the 
efficiency of that ideal cycle to approach the corresponding Carnot efficiency. 
The present techno-economic domination of practical Rankine-, Otto-, Diesel-, 
and Brayton-cycle engines, to the almost total exclusion of practical 
Ericsson-, Stirling-, or Carnot-cycle engines, is obvious. Yet the ideal 
thermodynamic efficiency of each of the three engine types that have thus far 
dominated industry is less than that of the Ericsson, Stirling, and Carnot 
engines. 

While the thermodynamic efficiency of any heat engine is of great 
interest, the work performed per cycle, per unit time, per unit volume, and 
per unit mass is of at least equal interest. The ideal efficiency of the 
Carnot cycle is the highest attainable at given temperatures for the source 
and the sink; however, the thermodynamic work accomplished per repetitive 
cycle of the ideal Carnot engine is relatively small. Thus, the efficiency of 
a Carnot engine as determined on the basis of shaft work produced can be 
severely degraded by even moderate mechanical inefficiencies. The thermo- 
dynamic work per cycle produced by the ideal Stirling engine is relatively 
much larger than that produced by the ideal Carnot engine; however, because 
the gas dynamic and mechanical friction losses in any practical Stirling 
engine are considerable, the shaft work per cycle produced by a practical 
Stirling engine must be appreciably lower than that predicted for the 
corresponding ideal Stirling engine. 

In order to obtain a high specific power, i.e., a high power per unit 
volume displacement, the thermodynamic processes employed by any heat engine 
must be carried out rapidly. This necessity virtually always leads to 
thermodynamic irreversibilities and a resultant reduction in engine 
efficiency. A high source temperature is desirable, not only because the 
Carnot efficiency is increased by increasing this temperature, but because the 
transport processes within the normally gaseous or vaporous thermodynamic 
medium are enhanced at higher temperature . 
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One reason for the widespread utilization of Otto- and Diesel-cycle 
engines has been that these are cyclic internal combustion engines in which 
the effective thermodynamic source temperature is quite high while the 
temporal mean temperature of the boundaries within the combustion chamber can 
be maintained at a level that can be tolerated by the materials constituting 
the boundaries. On the other hand, development of the steady-flow Brayton- 
cycle engine was delayed until fairly recently by the unavailability of 
materials that could withstand the source temperature required for reasonably 
efficient operation of the engine. 

Another reason for the predominance of Otto- and Diesel-cycle engines is 
that, because these are internal combustion engines, the heat added to the 
working medium is generated within the working space and does not have to be 
transferred through the solid boundary of the working space as is required in 
an external combustion engine. 

Because the Stirling engine is an external combustion engine, the 
following is true: 

(1) All of the heat entering the thermodynamic working space must be 
transferred through the heater tubes at a temperature at least 
slightly above the source temperature of the cycle. 

(2) That part of the boundary at the high-temperature end of the 
working space within the engine must operate continuously at 
approximately the cycle source temperature and at a mean pressure 
of 50 to 150 atmospheres. 

(3) All of the thermal energy transferred into the thermodynamic 
medium and not converted to work or lost by leakage through the 
boundaries must be rejected through the engine coolers. While 
this last requirement does not impose serious materials problems 
on the design of the Stirling engine, it does necessitate a heat 
rejection system possessing as much as three times the cooling 
capacity of the heat rejection system associated with an Otto- or 
Diesel-cycle engine of the same power rating operating between the 
same temperatures of source and sink. 
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SECTION III 


ATTRACTIVE CHARACTERISTICS OF THE STIRLING ENGINE 


The following characteristics of the Stirling-cycle engine make it an 
attractive candidate for solar thermal electric applications: 

(1) High thermal efficiency: The relatively high thermal efficiency 

of the Stirling engine is indicated by Figure 2-11. Even higher 
efficiencies appear possible if isothermal izers , low-drag 
regenerators, and hydrostatic bearings and seals can be success- 
fully incorporated into a carefully designed Stirling engine. 

(2) Dynamically balanced, valveless, closed-cycle operation: Many 

dynamically balanced Stirling engines, kinematic and free-piston, 
single and multiple cylinder, have been designed. There is 
scarcely a reason for considering a Stirling engine that is not 
dynamically balanced. Therefore, mechanical vibration of the 
engine assembly is minimal. 

Because the Stirling engine is valveless and, in the case of the 
solar-powered engine, has no intake from or exhaust to the ambient 
air, engine operation can be expected to be far quieter than Otto- 
or Diesel-engine operation. Even when a fossil-fuel combustor is 
employed as the heat source for a Stirling engine, the continuous 
nature of the combustion process can be expected to allow the 
engine to operate more quietly and with less pollution than a 
comparable Otto or Diesel engine. 

(3) Insensitivity of engine efficiency to engine size: Design 

analyses have indicated that little or no penalty in engine 
efficiency is incurred by Stirling engines of small size, e.g., 
one to twenty-five kilowatts output, as compared with large 
engines of up to two thousand kilowatts output (Ref. 11). 
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SECTION IV 


STIRLING ENGINE DEVELOPMENT PROBLEMS 


The following development problems have been identified for the Stirling 
engine : 

(1) Regenerator effectiveness: A regenerator accepts thermal energy 

from a working fluid during one phase of the Stirling cycle and 
returns the energy during another phase. 

Regenerator effectiveness is defined as, 



when E < 1, engine efficiency is reduced. 

In general terms, regenerator effectiveness increases and 
regenerator flow losses decrease with increasing regenerator void 
volume. Unfortunately, regenerator void volume is a part of the 
so-called dead volume of the engine, and increasing the dead 
volume reduces the pressure ratio and, hence, the performance of 
the engine. 

(2) Dead space: Dead space is that part of the working space that is 

not swept by the pistons. Stirling engine dead space includes: 

(a) Hot-end cylinder clearance volume 

(b) Heater internal volume 

(c) Regenerator void volume 

(d) Cooler internal volume 

(e) Cold-end cylinder clearance volume 

(f) Connecting duct volumes 

Foreseeable improvements in conventional Stirling engine design 
cannot reduce dead space below about one-half the minimum total 
working space volume. 

Dead space reduces work output per cycle, even with perfect 
regeneration . 

Based on first- and second-order analyses, dead space reduces 
engine efficiency even when regeneration is perfect. 
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(3) Expansion and compression processes are more nearly adiabatic than 

isothermal: Each of the primary volumes wherein the expansion and 

compression processes take place in the Stirling engine is 
enclosed by the cylinder end walls and the piston end. In 
conventional designs, the surface available for heat transfer to 
or from the enclosed gaseous medium is limited. Furthermore, at 
the usual cyclic rate of the engine, the distance from the 
enclosing surface into the gas to which the heating or cooling of 
the solid surface exerts any appreciable effect on the temperature 
of the medium is also severely limited. Thus, in conventional 
designs, the expansion and compression processes are much more 
nearly adiabatic than isothermal, and isothermal expansion and 
compression are required for ideal Stirling engine operation. 

In recent years, attempts have been made to design the expansion 
and compression spaces so as to make the thermodynamic processes 
occurring therein more nearly isothermal. The cylinder and piston 
ends have been fitted with mating sets of plane or axially 
symmetric interleaving fins comprising deep serrations or parallel 
plates. These so-called isothermalizers increase surface area and 
decrease the average distance of the medium from the bounding 
surface; thus, the concept rests on well established principles of 
heat transfer. It seems surprising that the concept has not been 
more thoroughly tested or more widely applied. This is even more 
surprising when one realizes that the application of isothermal- 
izers within the working space, coupled with condensing liquid 
metal heat transfer from the source at the hot end of the engine 
and of boiling heat transfer to the sink at the cold end of the 
engine, would allow reduction, and perhaps even elimination, of 
the dead volumes associated with the heater and cooler of the 
conventional design. Direct radiation of energy through a window 
into the hot end of the engine has also been suggested. 

(4) Stirling engine piston rings: Rings must operate in a hot, dry 

environment. The best currently used rings employ a teflon-filled 
inorganic matrix. Successful operation is measured in hundreds of 
hours. Higher-than-current operating source temperatures (>1350°F) 
will aggravate this problem. A gas bearing type of piston support 
and seal in the cylinder may be the best solution. However, two 
possible problems remain. 


(a) Low friction in the gas bearing implies a large clearance 
and large pumping losses, while small clearance in the gas 
bearing implies smaller pumping losses but higher piston 
friction losses. In contemporary gas bearings, the pumping 
losses predominate (Ref. 12). 

(b) The very close fit required for a satisfactory gas bearing 
complicates the problem of maintaining the correct clearance 
as the two mating parts heat and cool in transient operation 
and/or assume different temperatures in steady operation. 
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While the analogy is not perfect, it is sobering to consider the 
fact that, after more than fifty years of continuous development, 
properly lubricated, hard chrome rings in nitrided cylinder liners 
allow only about fifteen thousand hours of Diesel engine operation 
before overhaul or replacement is required. 

Piston rod seals: At least until very recently, the best Stirling 

engine rod seals have had an average operating life of a few 
hundred hours. This has been the most serious and stubborn 
operational problem of the modern kinematic Stirling engine . 

Rod seals must limit leakage of the working fluid from the work 
space to avoid frequent recharging of the engine. Tight sliding 
seals mean large friction losses . Rolling seals have exhibited a 
very inconsistent and, on the average, a very limited fatigue 
life . Rod seals must prevent any lubricant other than the working 
fluid from migrating into the work space because oil entering the 
working space is quickly pyrolyzed and deposited on the 
regenerators, severely reducing engine efficiency. When hydrogen 
is employed as a working fluid, some of any oil entering the 
working space reacts with the hydrogen to form methane, this 
admixing of methane with the hydrogen also reduces engine 
efficiency. 

As with the piston ring problem, a suitable gas bearing may 
provide a satisfactory solution — at least for the free-piston- 
type engine. 

Materials strengths at the high source temperatures and high mean 
pressure in the working volume: Even at the heater tube and 

heater head design temperatures that have been commonly employed 
in contemporary Stirling engines, viz., approximately 1350°F, 
failure of heater tubes and/or heater heads due to elevated- 
temperature creep or combined mechanical and thermal stresses 
occur occasionally. 

At the 1500 °f to 1600°F source design temperatures anticipated 
for the metallic engines now being developed and tested, these 
types of materials failures will continue to be a serious problem, 
even when the heater tubes and head are fabricated from expensive, 
predominantly non-ferrous, high-temperature alloys. 

At least one recently constructed prototype engine employed 
ceramic inserts in the heater head as a means for partially 
insulating the structural metallic head from the high-temperature 
medium. Currently, considerable interest is centered on the 
fabrication of complete heater tube and heater head assemblies 
from structural silicon carbide and/or silicon nitride. However, 
development of the techniques for constructing consistently sound 
subassemblies of these materials in a manner that will allow them 
to withstand high pressures, high temperatures, and large thermal 
transients has only begun. 


(7) Mechanical efficiency: The mechanical efficiency of the kinematic- 

type Stirling is seldom more than 85% and is usually 80% or less 

if a rhombic drive is employed. The mechanical efficiency of the 
free-piston-type Stirling is hardly more than 90% because of gas 
bearing and gas spring losses, and the relatively low linear 
alternator efficiency demonstrated to date ( < 93%) will offset the 
relatively higher mechanical efficiency of the free-piston engine 
unless linear alternator efficiency can be improved considerably 
in future designs. 

Friction losses in the piston and piston rod seals probably cannot 
be reduced very much unless hydrostatic sealing proves feasible. 
Similarly, piston pin friction and, in the case of the rhombic- 
drive kinematic engine, friction in the linkage pins that do not 
rotate continuously in one direction are probably irreducible 
unless hydrostatic bearings can be used in these locations. 

(8) Cooling system: The Stirling engine is an external combustion 

engine; thus, all heat entering the working space must be 
converted to work, removed by the cooling system, or leak to the 
surroundings. Of course, indiscriminate heat leakage must be 
limited in order to insure efficient operation, while the energy 
converted to work is limited by Carnot's principle; the remaining 
energy — roughly one-third to two-thirds of that supplied — must 
be removed by the cooling system. As a result, the required 
cooling system capacity for a typical current Stirling engine is 
roughly three times that of a current Otto or Diesel engine of 
equal power output. 

(9) The working medium: In order to provide high thermal capacity and 

high heat transport capability while, at the same time, incurring 
low turbulent flow losses, a low-density thermodynamic medium 
should be employed in an efficient Stirling engine. Hydrogen and 
helium are the two commonly used media. However, some concomitant 
problems are incurred with either of these choices. In any case, 
the mean working pressure in the usual Stirling engine is fifty to 
one hundred and fifty atmospheres, so an abrupt and potentially 
dangerous decompression is always possible except when the medium 
is withdrawn from the engine. 

Because the molecular weight of either hydrogen or helium is low, 
leakage of the medium from the working space is a problem unless 
hermetic sealing of the engine-generator assembly can be 
accomplished. Even then, hydrogen can diffuse through hot metal 
parts. This task appears now to be less formidable in the case of 
the free-piston Stirling engine than in the case of the kinematic 
Stirling engine. 

While the molecular weight of helium is approximately twice that 
of hydrogen, helium is monatomic while hydrogen is diatomic. 
Reference to Section II. B. 2., "Regenerative Stirling Cycle with 
Regenerator Dead Space," allows the conclusion that, unless 
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regeneration is perfect, the efficiency of this cycle will 
increase with increasing ratio of specific heats, k. Because 
k = 1.67 for a monatomic gas while k ~ 1.40 for a diatomic gas at 
the temperatures to be encountered in metallic Stirling engines, 
helium is to be preferred on this basis. On the other hand, given 
identical operating conditions, the flow losses can be expected to 
be higher when helium rather than hydrogen is employed as the 
medium. Two additional problems should be noted: Hydrogen is 

highly flammable, and helium, while non-flammable, is only 
available in limited quantities and at a higher cost. 

Consideration of some of the hazards noted above . strongly indicates 
that skilled mechanical technicians — equivalent in training and 
experience to licensed aircraft and engine mechanics — will be 
required for the satisfactory maintenance of Stirling engines. 
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SECTION V 


SUMMARY AND CONCLUSIONS 


The major conclusions of this investigation support and reinforce the 
already widely held opinion that the Stirling engine is strongly competitive 
both technically and economically as the power conversion subsystem of an 
advanced point-focusing distributed receiver solar thermal power system. 

In both magnitude and trend, brake efficiencies projected here are in 
satisfactory agreement with those estimated by United Stirling for advanced 
derivatives of their current P-40 and P-75 engines, e.g., a brake efficiency 
of forty percent (perhaps even one or two percent higher) at a heater head 
temperature of 1500°F and a cooler temperature of 140°F. Although the 
thermal efficiency of some current Stirling engines is quite high — superior 
to virtually any other practical engine operating in the same temperature 
range — - considerably higher efficiency is promised by such developing 
innovations in Stirling engine design as isothermalizers , high-heat-transfer , 
low-flow-loss heat exchangers, and hydrostatic seals and bearings. Minimiza- 
tion of working space dead volume is a key to high Stirling engine efficiency; 
therefore, it is imperative that the hot-end, regenerator, and cold-end dead 
volumes be minimized. Another key to high engine efficiency is the highest 
regenerator effectiveness consistent with the concomitant requirements of 
minimum regenerator void volume and low flow loss. Regenerator effectiveness 
becomes increasingly important as the ratio of heater head temperature to 
cooler temperature is increased; thus, higher heater head temperature implies 
more precise regenerator optimization if engine efficiency is to be increased. 

The Stirling engine appears to be well matched with the cavity receivers 
used in point-focusing distributed receivers in that highest combined 
efficiency appears to occur at no more than 1900°F, and this efficiency is 
only modestly reduced at receiver temperatures as low as 1600°F (Ref. 13). 

Of current major concern are the practical problems associated with assuring 
reliable, low-maintenance, long-life operation of practical Stirling engines 
that are cost-competitive with other engines, e.g., the gas turbine engine. 

The necessity for reducing or eliminating the use of expensive, high-temperature 
metal alloys that require strategically vulnerable alloying elements in Stirling 
engines is also clear but somewhat less pressing. If the durability, reli- 
ability, and serviceability of the kinematic Stirling engine is to be improved 
appreciably, the durability, reliability, and serviceability of piston rod 
seals, piston rings, and heater heads must be greatly improved. Therefore, it 
appears that, at least until considerable improvement is assured, advanced 
technological development activities should be concentrated on these components. 

While the promise that appears to be offered by properly designed advanced 
free-piston Stirling engines has been only outlined in this report, that promise 
should be the object of further investigation. 
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APPENDIX 


REFINED CALCULATION OF THE HEAT EXCHANGE IN REGENERATORS [0]* 

by H. Hausen 

(Translated by J. M. Bowyer) 


A theory of heat exchange in primarily ceramic regenerators is extended 
in such a way that even the rapid temperature changes immediately following 
the flow reversals and the deviations in temperature gradients at the ends of 
the regenerator can be considered. The equations and diagrams presented allow 
an exact calculation of regenerators, even under unusual conditions, and apply 
even in the limiting cases of arbitrarily thin or thick ceramics . 

In the language of the smelting industry, a regenerator is understood to 
be a reversible flow heat exchanger through which gas flows alternately in 
opposite directions and in which a mass of large thermal capacity periodically 
accepts and rejects the thermal energy that passes through it. 

In the case of ceramic regenerators operating at high temperature, e.g., 
the air preheater of a blast furnace, the heat exchange achieved depends, 
among other things, on the temperature variation within a cross section of the 
ceramic. On the basis of these temperature variations, Heiligenstaedt [1] and 
Rummel [2] have formulated equations for the calculation of a regenerator heat 
exchange coefficient or, respectively, a regenerator heat transfer 
coefficient. Heiligenstaedt [3] proceeded from the simplifying assumption 
that the gas temperatures do not vary with time; Rummel proceeded, in part, 
from empirical expressions. The author [4] has derived an equation that is 
based purely on theory and yet is simple and better fitted to the actual 
behavior (of a regenerator). However, as already emphasized in the (earlier) 
publication [4], calculation still neglects certain very rapid temperature 
changes which appear immediately after each flow reversal. Under the 
circumstances, even the deviating temperature gradient at the ends of the 
regenerator required a separate calculation for which a fully developed method 
was indicated at the conclusion of the report [5]. Under ordinary conditions 
the two effects mentioned (above) are usually small. However, in order to be 
able to consider them consistently correctly and yet quickly, an extended 
theory [6] is developed in what follows. In this way, it is possible to 
calculate very precisely the temperature gradient at all locations in a 
regenerator, as well as the heat transfer coefficient, even under unusual 
conditions . 


*Numbers in brackets represent Notes at the end of the Appendix. 
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Fundamental Concepts of the Previously Developed Theory 


The following investigations follow directly from the previously 
developed, simpler theory. Therefore, the basic concepts and the results of 
this (simpler) theory will be briefly discussed. We consider the temperature 
variation within a cross section of the ceramic that is perpendicular to the 
direction of flow of the gases (Figure A-l)*. To distinguish clearly, let the 
ceramic temperatures be denoted by t, the temperature of the warmer gas by 0, 
and that of the cooler gas by 0 ' . Because the outer surface of the ceramic is 
washed alternately by warmer and cooler gas, the ceramic temperature, t, 
oscillates constantly up and down between the two limits T max and T m £ n . 

These oscillations are larger near the outside surface than they are in the 
interior of the ceramic. How the oscillations under particular conditions 
proceed depends on the way in which the gas temperatures, 0 and 0', vary with 
time . 


On the basis of experimental results as well as on theoretical 
considerations concerning low-temperature regenerators, the simpler theory 
assumed that the gas temperatures at the cross section of interest varied 
linearly with time between each two successive flow reversals. With the help 
of the finite difference method of Binder and Schmidt [7], one can then obtain 
the approximate (spatial) temperature profiles at various fixed times as shown 
in Figure A-2. The figure shows the cooling of a 40-mm-thick ceramic during 
the cooling period, that is, during the time in which the cooler gas flows 
through the regenerator. The duration of this period, as well as (that) of 
the subsequent heating period, has been assumed to be 12.75 min. 

At the start of the cooling period, the temperature distribution, which 
was obtained from a separate consideration of the steady (oscillating) state, 
is represented by a parabola or parabolic curve that is concave upward. As a 
result of the initially very rapid cooling of the outer layers of the ceramic, 
this curve is deformed and, after about three minutes, is transformed into a 
parabola that is concave downward. Subsequently, the curve moves downward 
with unchanging form and constant speed. 

If one plots this temperature variation as a function of time and, 
indeed, for cooling and heating periods whose times are measured in opposite 
directions, one obtains Figure A-3. The dashed lines represent the gas 
temperatures, and the solid, closed line, the outside surface temperatures, 
t 0 and t 0 ' , of the ceramic during the heating and cooling periods. At 
other locations in the ceramic, the temperatures traverse closed lines but 
with smaller vertical amplitude. However, if one determines a local mean 
value for the temperature or T m ' of the ceramic across the entire width 
of the cross section by application of the approximate method at each instant, 
one is led to the remarkable result that the temporal variation of this mean 
value can be very well represented by a single straight line which is 


*Figures appear at the end of the Appendix following Notes. 
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traversed in opposite directions in the two periods (heating and cooling). 
This encompasses a further, more general perception. If the times of the two 
periods which, in the chosen representation, assume the same abscissa are 
called "mutually corresponding times," then for such times, corresponding to 
Figure A-3, 


T 


m 


m 


( 1 ) 


Figure A-3 shows furthermore that, aside from the very rapid temperature 
changes immediately after the flow reversal, even the outer surface 
temperatures t 0 and t q ' vary linearly during most of each period. 

Therefore, the assumption was made in the earlier theory that all temperatures 
in the ceramic varied with time in a completely linear way, as is illustrated 
in Figure A-4. On the basis of this assumption, a purely parabolic 
temperature distribution was obtained from the differential equation for heat 
conduction. Using Equation (1), this temperature distribution ultimately led 
to a relationship for the heat transfer coefficient for which K. Rummel [8] 
has introduced the following definition. Let T be the duration of the heating 
period and T' that of the cooling period. Then the heat transfer coefficient, 
k, is defined as the quantity of thermal energy accepted and then rejected by 
one of the regenerator heating surface in a complete period, T + T', 
divided by the temporal mean value of the temperature difference, 0 - 0' , 
between the warm and cool gases. The relationship for k derived on the basis 
of the completely linear variation of temperature with time is 


l 

k 


(T + T' ) 





( 2 ) 


where a and a 1 represent the thermal conductances during the warm and cool 
periods, and 6 and A represent the thickness and thermal conductivity of the 
ceramic. Accordingly, the quantity of heat transferred from the entire heat 
transfer surface, F, in one complete period, T + T', is 

Q = kF (T + T') (0 - 0 1 > M , (3) 


where (0 - 0') M represents the mean temperature difference between the two 
gases and k has the same value at all points within the regenerator. 

The Differential Equations for Heat Exchange 
in the Regenerator as the Basis of the Extended Theory 

The neglected (factors) of the previously developed theory that were 
mentioned in the introduction can now be more precisely identified. As shown 
in Figure A-3, not only at the outer surfaces, but also within the ceramic, 



very rapid temperature changes occur immediately following the flow reversal. 
However, sooner or later these variations tend toward the linear variation of 
temperature with time. These rapid temperature changes, hitherto neglected, 
cause a larger error the thicker the ceramic elements and the shorter the 
duration of the heating and cooling periods are. Furthermore, as a rule, the 
assumption of a completely linear distribution of the mean ceramic 
temperatures, x m and x m ' , as well as the relationship x ffl = x m ' for 
the mutually corresponding times of the two periods is very well justified in 
the large middle part of the regenerator. However, it is no longer valid at 
the ends of the regenerator where the variations of x m and x m ' are 
nonlinear and, as should be indicated more precisely, are bent, forming a 
hysteresis loop, vide Figure A-20. 

Now, in order to develop a theory that accounts for these two deviations 
numerically, we will consider all of the periodic temperature variations in 
the steady operation of the regenerator as induced temperature oscillations. 

We will analyze these oscillations in a manner similar to that (employed in 
analyzing) the fundamental and higher harmonic vibrations of a string. It 
will then become apparent that an exact calculation of the fundamental 
oscillation will enable formulation of the rapid temperature variations 
following the flow reversal. In contrast, it will only be possible to 
calculate the temperature deviation at the end of the regenerator by including 
the harmonics. The applicable equations for the fundamental and harmonic 
oscillations can be obtained from the differential equations governing heat 
exchange in the regenerator. Therefore, these differential equations should 
next be briefly discussed. 

The heat conduction equation is valid for the temperature changes inside 
a cross section of ceramic that is perpendicular to the direction of flow: 



where x represents the temperature of the ceramic at distance y from its 
left-side, outer surface (cf. Figure A-2), t the time, and a = A/(Yc) the 
thermal diffusivity. In turn, A represents the thermal conductivity of the 
ceramic (as noted earlier), Y its density, and c its specific heat. In 
Equation (4), the thermal conductance of the ceramic in the streamwise 
direction is neglected as unimportant [9]. 

Henceforth, let f be the heat transfer surface of the ceramic (elements) 
measured from the entrance point of the gas into the regenerator up to the 
cross section of interest, and let W be the thermal capacity of the gas 
flowing through the regenerator per unit time during the heating part of the 
cycle. As is easily determined, the following two differential equations for 
heat exchange in the regenerator are then valid: 
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and 



The meaning of these equations is immediately clear if one multiplies 
Equation (5) by -W and Equation (6) by Yc6/2, which represents the thermal 
capacity per unit area of the heated surface of the ceramic (considered to be 
a flat surface). 

Thus, Equation (5) expresses the fact that heat given up by the passing 
gas to unit surface area of the ceramic per unit time is equal to the heat 
transferred between the gas and the surface of the ceramic through the 
temperature difference 6 - x Q . Equation (6) states that the same quantity 
of heat is transferred from the gas to the regenerator (represented by the 
right member of Equation (5) multiplied by W) and, in consequence, the mean 
temperature of the ceramic, x m , increases by (9x m /9t)f P er un i t time. 

Equations (5) and (6) also apply to the cooling periods if 0, x Q , x m , a, 
and W are replaced by 0 ' , x 0 ', x m ' , a', and W' and if heat transferred 
from the ceramic surface to the gas is now counted as positive. 

In order to find solutions of these equations which describe the steady 
state of the regenerator, the further condition must be fulfilled. At the end 
of a complete cycle consisting of a heating and a cooling period, the 
temperature at every point in the ceramic must be the same as it was at the 
start of that cycle. This condition is called the flow reversal condition. 
There are infinitely many solutions which satisfy this condition. These are the 
eigenfunctions, and, at times, are indicated by an "eigenvalue," k, which may 
be any integer of the infinite series that extends from 0 to », 

As will be shown, the eigenfunction for k = 0 represents the fundamental 
oscillation, while the higher eigenfunctions k > 1 yield the harmonics. 

For simplicity, only the eigenfunctions for the case that the heat 
capacities of the two gases streaming through the regenerator are temperature- 
independent and the durations of the heating and cooling periods are equal 
will be discussed. If, as earlier, W and W' represent the heat capacities of 
the gases per unit temperature, and T and T' , the duration of the heating and 
cooling periods, this (the foregoing limitation on eigenfunctions) means we 
will limit ourselves to the case WT = W'T'. However, it will be briefly 
indicated in conclusion that the results of these considerations remain 
applicable and sufficiently accurate, even in general cases. 


The Fundamental Oscillation (Zeroth Eigenvalue) 

In most earlier works regarding the calculation of ceramic regenerators, 
it was assumed as obvious that the temperatures in the lengthwise direction of 
a regenerator varied just as they do in a recuperator, i.e., as in steady flow 
through a tubular heat exchanger. In this analogy, the temperature distribu- 
tion in the regenerator is considered at a definite time or at a time corres- 
ponding to the middle of a (heating or cooling) period. Among the known 
solutions of the differential equations, one can always be found that 
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represents this simplest longitudinal temperature distribution, and it can be 
demonstrated that this solution is obviously the zeroth eigenfunction (k = 0) 
that corresponds to the fundamental oscillation [10]. Thus, the zeroth 
eigenfunction can also be defined as that which yields the same longitudinal 
temperature distribution (in the regenerator) as is found in a recuperator. 
Then if, hereafter, we seek the zeroth eigenfunction for (the case) WT = W'T', 
we must, for the analogy, refer to a recuperator in which the thermal 
capacities of the two gases are likewise equal. It is known that, in this 
case, a completely linear temperature distribution is established in the 
recuperator as long as W, W', a, and a’ are constant, cf. Figure A-5. 
Therefore, the zeroth eigenfunction must present, at every instant, a linear 
temperature distribution in the longitudinal direction of the recuperator. 
Because the longitudinal coordinate was chosen to be the heat transfer 
surface, f, extending up to the cross section of interest, hereafter 



= constant 


(7) 


will be valid for the gas temperature, 0, as well as for the temperature of 
the ceramic, t, at any particular time, t. 

Equation (7) immediately allows important conclusions to be drawn from 
the differential equations if, beside a and W, the parameters A, Y, c, and <5 
are also assumed constant. First, from Equation (5), it follows that 
t 0 - 0 “ constant, that is to say, the temperature difference between the 
outer surface of the ceramic and the gas remains constant during the entire 
(heating or cooling) period. Furthermore, from Equation (6) we obtain 


= constant. (8) 

f 

From the linear variation (of x m ) in the longitudinal direction [and 
Equation (8)], it follows that the mean temperature of the ceramic at each 
cross section varies linearly with time. This rule that was viewed earlier as 
only a very good approximation is thus proved to be an exact requirement of 
the zeroth eigenfunction when WT = W'T 1 . Consequently, for mutually 
corresponding times of the two periods, Equation (1), 



is exactly satisfied by the zeroth eigenfunction. Recognizing this, the mean 
temperatures of the ceramic, x m and x m ', become of fundamental importance 
to further considerations. Certainly, they form a simple framework upon which 
the essentially complicated behavior of the individual temperatures of the 
ceramic, x and x', and the gas temperatures, 0 and 0' will be constructed. 



The Temperature Behavior According to the Zeroth Eigenfunction 


Figures A-2 and A-3 already provide an approximate demonstration as to 
how the gas temperatures and ceramic temperatures behave when restricted to 
the zeroth eigenfunction. In order to be able to calculate exactly the 
temperature behavior under widely different conditions, we employ the 
following solution of the differential equation [Equation (4)]: 


T = C + 


n — 1 


-6 at 


n 


cos g 


n (* “ I) 


(9) 


where C, Bq, and g n represent initially arbitrarily chosen constants and 
(9T m /9t)f is invariant in consequence of Equation (8). The lower case 
letter n represents a positive integer. The first term following C expresses 
the fact that x m should vary linearly with time. The second term presents 
the parabolic temperature distribution as it tends to develop toward the end 
of the period, cf. Figure A-2, A-7, or A-18. The last term, involving the 
cosine, yields the deviations from the parabolic distribution which are 
important immediately after the flow direction is reversed. Incidentally, 
because this member is a function of y, it is inferred that the solution must 
be symmetric with respect to the center of the ceramic, y = 6/2. 

In order to maintain the linear behavior of x m , we determine the 
values of B n so that the mean value evaluated over the entire thickness of 
the ceramic, 6, causes every cosine term to disappear; thus, these members 
make no contribution to x m . Therefore, we set 


J. 

6 



cos 



and thereby obtain 



Similarly, the mean value of the parabolic term disappears if 6^/6 is 
subtracted from y(6 - y). Because it is still possible arbitrarily to choose 
the sign of B n by considering Equations (6) and (10), Equation (9) can also 
be written as 


1 = T » ‘ ft <e ' V L y - y > - r. 


00 / 2 nT l ^ at 

+ B n e * cos (^niT (11) 

n=l 


We associate this equation with the heating period; correspondingly, for the 
cooling period we obtain 


T = T 


m 


- Xs (e * - 


x ') 
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_y - y> - 1-] 



(12) 


We can further transform the last equation if we consider that, per unit 
surface area, the same quantity of heat, q, is exchanged in the cooling period 
as is exchange in the heating period. Thus, 


q = a' (t ' - 6') T' = a (8 - t ) T 
o o 


(13) 


applies, and, thus, Equation (12) is transformed into 
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In each period, the time t or t ' should be redefined as zero at the instant 
flow reversal begins. 

The values of Bq and Bjj' are obtained from the flow reversal 
condition whereby the ceramic temperature, x, at the beginning of the heating 
period should have the same value as at the end of the cooling period, x', and 
vice versa. If we determine the temperature at the start of the heating 
period (t = 0) from Equation (11) and the temperature at the end of the 
cooling period ( t ' = T ' ) from Equation (14), then equating (these two 
equations) and recalling Equation (1), 



On the left side of this equation is a Fourier series with the yet-to-be- 
determined coefficients, 


B - B ' e ' 6 
n n 




If the right side of Equation (15) is defined as f(y), then according to the 
known formula for determining the Fourier coefficients, 


B - B ' e 
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and, performing the integration, 
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A similar consideration of conditions at the end of the heating period 
and the beginning of the cooling period yields 
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If Equations (16) and (17) are solved for B n and B n ' and then B n 
is substituted into Equation (11), the most general form of the zeroth 
eigenfunction for WT - WT' is 



In the above equation, by comparison with Equation (6), x m is given by 



(19) 


where (x m ) a represents the value of T m at the start of the period. 

Using Equations (18) and (19), calculate the complete temperature behavior for 
the period of interest and, by a corresponding application (of the method) 
also calculate the complete temperature behavior in the immediately following 
period. The gas temperature 0 is obtained from Equation (18) by calculating 
the ceramic surface temperature (y = 0) and adding it to the temporally 
invariant quantity 0 - t q . 


Comparison of Theory with Experiment 

Although the accuracy of the relationships found follows from the 
rigorous, exact derivation carried out (above), it furthermore may be shown 
that it is also confirmed by experiment. The temperature behavior calculated 
by Equations (18) and (19) in an 80-mm-thick ceramic, for which the period 
durations T = T' = 1 h are assumed and for WT = W'T' , is shown in Figures A-6 
and A-7 . The curves in Figure A-6 give the temporal temperature distribution 
at a specified cross section of the ceramic at its outer surface (y = 0), at 
distances of 15 and 25 mm from the outer surface, and at its mid-plane 
(y = 40 mm). (See also Figure A-3.) The linear behavior of the mean 
temperatures of the ceramic, t m and T m ' , is indicated by the dot-dash 
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line. As in Figure A-2, Figure A-7 shows the temperature distribution in the 
ceramic's cross section at different times, t. For comparison with the curves 
in Figure A-6, the measurements of Schumacher [13] are presented in 
Figure A-8. The experimental regenerator has the same characteristics as 
those relating to Figure A-6. The agreement between calculation and 
experiment is, to the accuracy of the experiments, very good. In particular, 
it can be recognized that even the rapid temperature changes immediately 
following the flow reversal are correctly represented by the newly developed 
equations. The good agreement also extends to the temporal behavior of the 
gas temperatures, 6 and 6'; calculated results are shown in Figures A-9 
through A-12, and the experimental results of Schumacher are shown in Figures 
A-13 through A-16. Also, in contrast to the assumptions imposed in Figure 
A-3, the temperature trend (presented) here is curved immediately after the 
flow reversal and becomes monotonically more nearly linear toward the end of 
the period. When comparing results, note that the measurements across an 
80-mm-thick ceramic are presented in Figures A-13 through A-16 while the 
calculated curve in Figure A-12 corresponds to a 200-mm-thick ceramic. 

(Compare also Figures A- 17 and A- 18.) 

How the temperature distribution in the ceramic varies if the ceramic is 
thicker or the period shorter should now be further explained. The thicker 
the ceramic, the longer the time (required) for transition from parabolic or 
quasi-parabolic curve concave upward to the concave downward form. For very 
thick ceramic this (transition) can require so much time that, even at the end 
of the period, the fully parabolic form is not completely achieved. 
Furthermore, the deviation from this form becomes larger as the period 
shortens. Figures A-17 and A-18 show the result of a precise calculation for 
a 200-mm-thick ceramic under otherwise identical conditions as (those imposed) 
in Figures A-6 and A-7 and where, in particular, the duration of a period is 
assumed to be 1 h. One recognizes that the middle (region) of the ceramic 
participates very little in the temperature oscillations and, for this reason, 
in the heat exchange. Assuming a constant time for the period, the 
temperature oscillations in the middle (region) of the ceramic are virtually 
non-existent beyond a certain thickness of the ceramic. A further increase in 
the thickness of the ceramic can then no longer affect the heat exchange. 
Because of this, the heat transfer coefficient must also be independent of 
ceramic thickness. 


The Heat Transfer Coefficient, k„, Corresponding to the 
Zeroth Eigenfunction 

Once the temperature behavior is determined from Equations (18) and 
(19), the heat transfer coefficient, k Q , that corresponds to the zeroth 
eigenfunction can be determined without special difficulty. 

First, the temporal mean value of the ceramic's external surface 
temperature during the heating period, (t 0 ) m , is determined by setting 
y = 0 in Equation (18) and integrating (the result) with respect to time, t, 
from 0 to T. Thus is obtained 


(t ) - (t ) = ^ (e - t )<j) , 

o M m M 6A o 7 


( 20 ) 
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where (T m ){^ represents the temporal mean of T m and is an abbreviation 
for 


<t> = 1 




( 21 ) 


In a similar way, for the cooling period, 


(t') m - 
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( 22 ) 


is obtained. 

Now the temporal mean value, (6 - 0 ')m> representing the difference of 
the gas temperatures, can be calculated with Equations (20) and (22). Because 
0 - t q = constant, x 0 ' - 0' = constant, and (T m ) M = (T m ' ) M [from Equation (1)], 
(0 - 9 ' ){4 can be determined from 

(T ') - (T ') 1 + (T ' - 6') . 

m M o mj o 

(23) 



By substituting Equations (20) and (22) into Equation (23), recalling 
Equation (13), 




is obtained. 

From the definition of the heat transfer coefficient, for each m^ of 
heating surface over a complete (regenerator) cycle, T + T' , the quantity of 
heat, q [= Q/F, cf. Equation (3)], absorbed and again given up by the ceramic 
is 


q=a(0-x)T=k (0- 0')^ (T + T'). (25) 
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From this (last equation), employing (0 - 0 ' )m from Equation (24), the 
following formula for the heat transfer coefficient, k Q , is obtained: 


l 

k 


(T + T' ) 



( 26 ) 


This equation differs from the previously derived formula [Equation (2)] only 
by the appearance of the factor <t>. This factor expresses the effect of the 
very rapid temperature changes following flow reversal. While, as a result of 
neglecting this effect, the earlier equation was the more inaccurate the 
thicker the ceramic and the shorter the duration of the complete period, 
Equation (26) remains exact in the realm of solutions expressible by the 
zeroth eigenfunction for infinitely thick ceramic and for arbitrarily long 
durations of the period. 

The appearance of the function <t> [in Equation (26)], which can be 
determined from Equation (21), appears at first to make the calculation of 
k Q unnecessarily difficult. However, it is possible to replace Equation 
(21) by the following very accurate approximate formulae: 


£or h (f + ?■) s 10 - ♦ 

f ° r I ; (? * ?-) > i°, 4 

It is even simpler to obtain <j> from 
is presented graphically as a function of 6 2 /2T or, as the case may be, 

(<5 2 /2a) (1/T + 1/T'). The solid curve is valid for T = T', while the dashed 
curves are valid for T = 2T' or T = 1/2 T' or for T - 4T' or T = 1/4 T'. The 
dashed lines show that, even for T = 4T 1 or T = 1/4 T' , the deviations from 
the values for T = T' are very small. Therefore, in almost all practical 
cases, it will be possible to determine <p with sufficient accuracy by means of 
the graph (of Figure A-19) or from Equations (27) and (28), which correspond 
very closely to the curves. 
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Figure A-19 wherein d> from Equation (21) 


The Limiting Cases of Very Thin and Very Thick Ceramics 

For very thin ceramics, 6 2 /(aT) or (6 2 / 2a) (1/T + 1/T') is very 
nearly 0; then, according to Figure A-19, 4> will be equal to one, and Equation 
(26) will be transformed into Equation (2). Thus, in terms of the zeroth 
eigenfunction solution, the thinner the ceramics, the more accurately 
applicable this simpler relationship becomes. Furthermore, Equation (28), as 
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well as an exactly determined limiting value (viz. appendix), shows that 
Equation (21) tends toward the expression, 


4 > 


2.142 HI a 



(for very thick ceramics), 


( 29 ) 


as the thickness, 6, of the ceramic increases. Correspondingly, Equation (26) 
assumes the following form, 


i_ = ( T + T') 
o 
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(for very thick ceramics). 


(30) 


The thickness of the ceramics no longer appears in this expression. Thus, the 
already mentioned requirement that, for very thick ceramics k Q must be 
independent of ceramic thickness, is fulfilled. 

The asymptotic representation [Equation (29)] is so good, that from 
(S 2 /2a) ( 1/T + 1/T' ) = 10 on up, it is almost always possible to apply 
Equation (30) (with) satisfactory accuracy. 


Influence of the Dust Layer 


The influence of the dust layer, which was first evaluated by Schack 
[14], can be calculated very simply through a slight extension of the theory 
that has been developed to this point. 

A dust layer of thickness 6 Q and thermal conductivity A Q is to be 
found on each of the two sides of the ceramic. If 6 0 is so small that the 
thermal capacity of the dust layer can be neglected by comparison with the 
thermal capacity of the ceramic, then the entire heat given up by the gas must 
be transferred through the dust layer. Therefore, the heat arriving at the 
ceramic outer surface from the gas has not only overcome a thermal resistance 
proportional to 1/a but also the resistance of the dust layer, which is 
proportional to & Q /\ Q . Thus, it is necessary to replace 1/a by 1/a + 

6 0 /X 0 and 1/a' by 1/a' + tS 0 A 0 in the earlier equations. In this way, 

Equation (26) for the heat transfer coef f icient, k Q , is transformed into 


1— = (T + T') 
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o 



(31) 
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By comparison with Equation (26), the essential factor expressing the 
influence of the dust layer in Equation (31) is, accordingly, the added term 
6 0 /A 0 . In this respect it is to be borne in mind that the heat transfer 
coefficients, a and a' can have somewhat different values in the case of heat 
transfer to the dust layer than (exists) in the case of transfer to the 
dust-free outer surface of the ceramic. 


The Temperature Behavior at the Ends of the Regenerator 
and Its Influence on the Heat Exchange (Process) 

As a rule, the zeroth eigenfunction very accurately represents the 
temperature changes within the inner regions of a regenerator. As already 
indicated, the processes (occurring) at the ends of the regenerator deviate 
from the trend of the zeroth eigenfunction. In accordance with the usual 
practical case, we assume that, as the case may be, the warmer as well as the 
cooler gas enters the regenerator at constant temperature. Then, for example, 
at the warm end of the regenerator, the temperature behaves basically as shown 
in Figure A-20. Again, the times are drawn in opposite directions in the 
warming and cooling periods. The uppermost curve, 0 = 6 , represents the 
temporally invariant entrance temperature of the warmer gas. The upper solid 
curves show the behavior of the outside surface and the mean temperature of 
the ceramic in the heating period. During the cooling period, the 
corresponding temperatures are 0', t q ' and T m ' . The temperature of the 
cooler gas, 0', that exits here varies less with time than it does in the 
inner region of the regenerator. It is most noteworthy that the curves for 
ijp and T m ' are not coincident as in Figure 3, for example. The zeroth 
eigenfunction is unable to represent either the constant entrance temperature 
of the gas or the curved trend of i m and x m ' . To obtain an exact 
calculation, it is necessary to introduce the higher (order) eigenfunctions. 

Certainly, it is possible to perform the calculation without the higher 
eigenfunctions if one ascertains the temperature behavior with respect to 
deviations from the zeroth eigenfunction by a temporal or other approximation 
process [15]. The effectiveness of the regenerator, which will be discussed 
later, can be calculated in this way. For practical numerical calculations, 
such methods may even be preferable. However, the physical relationship 
becomes clearest when one considers the higher eigenfunctions and investigates 
their influence on the heat exchange (process). 


On the Heat Transfer Coefficient Based on the 
Central Temperature of the Ceramic 

Because the higher eigenfunctions behave in an essentially more complex 

way than the zeroth eigenfunction and because at most they play only the role 

of correction terms, it is possible to avoid them in a complete and exact 

determination of the local temperature variations in a ceramic cross section 

at each instant. Furthermore, it suffices to assume the difference between 

T o and T m > and therefore between 0 and i m , at every instant to be equal 

to the time-averaged heat transferred per unit time calculated by means of 

the zeroth eigenfunction [16]. The temporal mean value of t. - i 

o m ° 
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determined from Equation (20). If the quantity 0 - T 0 , which is invariant 
in the case of the zeroth eigenfunction, is added to both sides of this 
equation, then 


(9 - 



. a . (6 - t ) 
o 


is obtained as the temporal mean value of 9 - T m . Furthermore, because 
q = a *(0 - x 0 )-T represents the quantity of heat transferred per unit 
area during one period, it is also possible to write, in place of the 
foregoing equation, 


q = a (6 




(32) 


if, for abbreviation, we set 



(33) 


According to Equation (32), we can consider a as a heat transfer coefficient 
based on the (temporal mean) temperature at the middle of the ceramic (cross 
section). 

With this heat transfer coefficient, a, we will now calculate (dq) in 
the case of the higher eigenvalues at each point in time by writing 


dq = ot(0-x )dt, 

n m ’ 


(34) 


for the quantity of heat transferred during the short time dt in place of 
Equation (32). 

The error that is made by applying Equations (33) and (34) to the higher 
eigenfunctions plays only a very small role in the end results as the 
following consideration shows. It is unthinkable that, through the value of 4> 
in Equation (33), the influence of the very rapid temperature changes 
immediately after the flow reversal will be, to some degree, uniformly 
distributed through the total duration of the period, particularly for the 
calculation of the heat transfer coefficient (which is our ultimate goal). 
From this point of view, it appears at first that very noticeable errors 
(will) occur in the case of the higher eigenfunctions because the gas and 
ceramic temperatures as determined by these functions behave essentially 
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differently than the zeroth eigenfunction. Therefore, it cannot be expected 
that the higher eigenfunctions will be well approximated by Equation (34). 
However, for the accuracy of the final result, only the magnitude of the error 
in suitably formed sums, yet to be more precisely discussed, of the higher 
eigenfunctions is significant. It was shown in an earlier publication [17] 
for the case <j> = 1 that this error is small and, even in unfavorable cases, 
affects the total quantity of heat transferred in the regenerator by less than 
1%. In the case <P < 1, even smaller errors are to be expected because then, 
according to Equation (20), even the temperature difference t 0 - x m is 
smaller, and an inaccuracy in calculating it is correspondingly less important. 

The introduction of ot in the case of the higher eigenfunctions makes it 
possible to calculate (the heat transfer) as if the temperature of the 
regenerator core at a given cross section and given instant had the uniform 
temperature, x m , and as if the temperature difference, 0 - i m , 
corresponding to the heat transfer to the outer surface, were at one's 
disposal. This limiting case of uniform temperature over the entire cross 
section is applied with very good accuracy to the regenerators of 
low-temperature technology because these regenerator cores are fabricated from 
thin plates of high thermal conductivity instead of ceramic. In his first 
report on regenerators [18] the author derived the equations for all the 
eigenfunctions corresponding to this case. The calculation with a and x m , 
therefore, allows these equations also to be used without change in current 
considerations. Even if the already mentioned approximate_methods are used 
instead of the higher eigenfunctions, the introduction of a is necessary 
because even these processes have been developed for uniform temperatures in 
the ceramic cross section. 

It should also be mentioned that, in the case of the dust-covered 
ceramic outer surface, Equation (33), corresponding to the foregoing 
considerations, is to be replaced by 


1 


a 


1 

a 



(33a) 


The Higher Eigenfunctions 


The paper continues, but translation is terminated at this point. 
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NOTES 


0. Extended edition of a contribution of the author to the Convention of the 
VDI - Committee for Thermal Research on 13 October 1941 in Dresden. See 
also the contribution to the convention by E. Schmidt in Z. VDI , 86 , 9/10 
(1942) pp. 135-138. 

1. W. Heiligenstaedt , "Die Berechnung von Warmespeichern, " Arch ■ 
Eisenhiittenw. , 2^ (1928/29) pp. 217-222; also W. Heiligenstaedt 
''Regenoratoren , Rekuperatoren , Winderhitzer , " Leipzig (1931). 

2. K. Rumrael, "Die Berechnung der WSrmespeicher auf Grund der 
Warmedurchgangzahl , " Stahl und Eisen , 48 (1928) pp. 1712-1715. 

3. In an earlier work, Heiligenstaedt has also discussed the case of the 
temporal variation of the gas temperatures, cf. W. Heiligenstaedt, 
"Beitrag zur Berechnung von Winderhitzern und Regeneratoren , " Piss . 

Aachen (1924) or p. 169 of the book mentioned in Note 1, above: 
"Regeneratoren, Rekuperatoren, Winderhitzer." 

4. H. Hausen, "Berechnung der Steintemperatur in Winderhitzern," Arch . 
Eisenhiittenw. , 12 (1938/39) pp. 473-480. 

5. H. Hausen, op. cit. p. 479 (Abschnitt: Erweiterung des 
Berechnungsverfahren) . 

6. I thank Dr. A. Schack in Diisseldorf for the encouragement to develop the 
extended theory presented here. Dr. Schack mentioned Equation (18) 
(corresponding to Equation (2) of this discussion) derived in my 
publication noted in footnote 4, in his book Der industrielle 
Warmeaustausch 2nd ed. Diisseldorf (1940) pp. 224 and 245. This has led 
to a long correspondence with me concerning questions raised there. As 
it has turned out, his objection to my equation is only that it does not 
apply even to thin ceramics in those cases where the deviating 
temperature behavior at the ends of the regenerator exerts a strong 
influence but is not considered in the calculation. 

7. See for example the publication of E. Schmidt in "Beitrage zur 
technischen Mechanik und technischen Physik," Festschrift August Foppl 
zum 70 Geburtstag Berlin (1924) or H. Grober and S. Erk, Die Grundgesetze 
der Warmeubertragung , 2nd ed . , Berlin (1933) p. 92. 

8. K. Rummel, op. cit.; K. Rummel and A. Schack, Stahl und Eisen , 49 (1929) 
pp. 1300-1311; Arch. Eisenhiittenw. , 2^ (1928/29) pp. 473-479; see also 

A. Schack, Arch. Eisenhiittenw, , (1928/29) pp. 481-486 and K. Rummel, 
Arch. Eisenhiittenw. , 4 (1930/31) pp. 367-374. 

9. If desired, this influence can be subsequently and separately 
considered. The author has presented a simple method suited to these 
conditions at the Convention of the Study Group "Kaltentechnik" of the 
VDI in Prague on 23 May 1941. See also K. Nesselmann, "Der Einfluss der 
Warmever luste auf Doppelrohrwarmeaustauscher ," Z. ges, Kalte-Ind., 35 
(1928) p. 62. 
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10 . 


For the simplest case, W = W* and T = T' ; this follows from Equations 
(46) and (47) of the following publication: H. Hausen, Uber die Theorie 

des Warmeaustausches in Regeneratoren, " Z.a ,M.M. , j) (1929) pp. 173-200. 
Compare Equations (35) and (40) of the above article. 

11. Mathematically, it is worth mentioning that the two members of Equation 
(9) that include (9T m /3t)f can also be considered to be the zeroth 
member of the cosine series, as can be shown easily through a limiting 
process by allowing g n — *" 0 as ^ — ► » . 

12. viz., e.g., B. W. Hort and A. Thoma, Die Dif ferentialgleichungen der 
Technik u. Physik , A. Barth, Leipzig (1939) p. 342. 

13. K. Schumacher, "Grossversuche an einer zu Studienzwecken gebauten 
Regenerativekammer , " Arch. Eisenhiittenw. , 4 (1930/31) pp. 63-74. 

14. A. Schack, "Uber den Einfluss des Stabbelages auf den Wirkungsgrad von 
Regeneratoren," Z. Techn. Fhys. , 9^ (1928) pp. 390-398; Arch . 
Eisenhiittenw, , 2_ (1928/29) pp . 287-292. 

15. H. Hausen "Naherungsverfahren zur Berechnung des Warmeaustauches in 
Regeneratoren," Z.a.M.M. , 11 (1931) pp. 105-114; see also H. Hausen, 
"Warmeaustauch in Regeneratoren," Z. VDI , 73 (1929) p. 431. 

16. The following consideration was already (though) less precisely and less 
generally presented in the earlier publication (cf. Notes 4 and 5). 

17 H. Hausen, op. cit. (cf. Notes 4 and 5) Figure 16. 

18. H. Hausen, "Uber die Theorie des Warmeaustauches in Regeneratoren," 
Z.a.M.M. , 9 (1929) pp. 173-200. 
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NOTES: 

e~ TEMPERATURE OF THE 
WARMER GAS 

Q'~ TEMPERATURE OF THE 
COOLER GAS 



Figure A-l. Oscillation of the Ceramic Temperature of a 
Regenerator Between T max and r m ^ n 



Figure A-2. Temperature Behavior in a Ceramic of Thickness 
6 During the Cooling Period 
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LEGEND: 

HEATING PERIOD: 

T o - OUTER SURFACE TEMPERATURE 
OF THE CERAMIC 

T m ~ LOCAL MEAN OF THE 
CERAMIC TEMPERATURE 

t ~ TIME 

COOLING PERIOD: 

~ OUTER SURFACE TEMPERATURE 
OF THE CERAMIC 

r m ~ LOCAL MEAN OF THE CERAMIC 
TEMPERATURE 

t' ~ TIME 


Figure A-3. Actual Temperature Behavior in the Ceramic 


LEGEND: 

HEATING PERIOD: 
t g - OUTER SURFACE TEMPERATURE 
OF THE CERAMIC 

~ LOCAL MEAN OF THE 
CERAMIC TEMPERATURE 

t - TIME 

COOLING PERIOD: 
r o ~ OUTER SURFACE TEMPERATURE 
OF THE CERAMIC 

Tm ~ LOCAL MEAN OF THE CERAMIC 
TEMPERATURE 

t'-TIME 


Figure A-4. Previously Assumed Temporal Temperature 
Behavior in the Ceramic 
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Figure A-5. Temperature Behavior in a Recuperator for Equal Thermal 
Capacity of the Two Gases (W = W') 
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Figure A-6. Calculated Temperature Behavior in an 80-mm-Thick Regenerator 
Ceramic (Cross Section) (See also Figure A-7) 
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NOTES: 

DURATION OF COOLING PERIOD: 1 h 

y~ DISTANCE FROM ONE OF THE 
TWO CERAMIC OUTER SURFACES 



y, mm 

Figure A-7 . Calculated Temperature Behavior of an 80-mm-Thick Regenerator 
Ceramic (Cross Section) (See also Figure A-6) 



Figure A-8. Measured Temperature Behavior in an 80-mm-Thick Regenerator 
Ceramic (According to K. Schumacher) 
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LEGEND: 

3~ CERAMIC THICKNESS 

r~ DURATION OF THE 
HEATING PERIOD 

t'~ DURATION OF THE 
COOLING PERIOD 





Figures A-9 through A-12. Calculated Temporal Behavior of the Gas Temperatures 

in a Regenerator. Heat transfer coefficients in 
the heating and cooling periods are 
a = a' =20 kcal/(m 2 h°C) . 
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LEGEND: 

8~ CERAMIC THICKNESS 

t~ DURATION OF THE 
HEATING PERIOD 

t'~ DURATION OF THE 
COOLING PERIOD 


Figures A-13 


TEMPERATURE 




Measured Temporal Behavior of the Gas 
Temperatures in a Regenerator 
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NOTES: 

DURATION OF 
HEATING AND COOLING 
PERIODS: 1 h EACH 

y~ DISTANCE FROM THE 
CERAMIC OUTER SURFACE 


Figure A-17. Temperature Behavior in a 200-mm-Thick Regenerator 
Ceramic (See also Figure A-18) 


NOTES: 

DURATION OF 
COOLING PERIOD: 

1 h 

/ ~ DISTANCE FROM 
ONE OF THE TWO 
CERAMIC OUTER 
SURFACES 



y, mm 

Figure A-18. Temperature Behavior in a 200-mm-Thick 

Regenerator Ceramic (See also Figure A-19) 
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Figure A-19. <j> as a Function of (S^/aT or (g2 /2a) (l/T + 1/T') as the 

Case May Be. (Required for the calculation of the heat 
transfer coefficient, k Q , for regenerators using 
Equation (26).) The value of Equation (21) can be 
obtained from the figure. 

LEGEND: 

8 ~ TEMPERATURE OF THE GAS 
IN THE WARMING PERIOD 

^'-TEMPERATURE OF THE GAS 
IN THE COOLING PERIOD 

t 0 ~ OUTER SURFACE TEMPERATURE 

OF THE CERAMIC IN THE 4 

HEATING PERIOD l 

LU 

-OUTER SURFACE TEMPERATURE 3 

OF THE CERAMIC IN THE COOL- < 

ING PERIOD £ 

r m ~MEAN TEMPERATURE OF THE | 

CERAMIC (OVER A CROSS 
SECTION) IN THE HEATING 
PERIOD 

— MEAN TEMPERATURE OF THE 
CERAMIC (OVER A CROSS 
SECTION) IN THE COOLING 
PERIOD 

Figure A-20. Temporal Temperature Behavior at the Warm End 
of the Regenerator 
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